Convection in Astrophysics 

Proceedings IA U Symposium No. 239, 2007 

F. Kupka, I. W. Roxburgh & K.L. Chan eds. 



© 2007 International Astronomical Union 
DOI: 00.0000/X000000000000000X 



Effects of Rotation and Input Energy Flux 
on Convective Overshooting 

Petri J. Kapyla 1 2 , M.J. Korpi 1 , M. Stix 2 
and I. Tuominen 1 

Observatory, P.O. Box 14, FI-00014 University of Helsinki, Finland 

2 Kiepenheuer-Institut fur Sonnenphysik, Schoneckstrasse 6, D-79104 Freiburg, Germany 

email : pkapy la@kis .uni-freiburg.de 

Abstract. We study convective overshooting by means of local 3D convection calculations. 
Using a mixing length model of the solar convection zone (CZ) as a guide, we determine the 
Coriolis number (Co), which is the inverse of the Rossby number, to be of the order of ten or 
larger at the base of the solar CZ. Therefore we perform convection calculations in the range 
Co = 0...10 and interpret the value of Co realised in the calculation to represent a depth in 
the solar CZ. In order to study the dependence on rotation, we compute the mixing length 
parameters «t and a u relating the temperature and velocity fluctuations, respectively, to the 
mean thermal stratification. We find that the mixing length parameters for the rapid rotation 
case, corresponding to the base of the solar CZ, are 3-5 times smaller than in the nonrotating 
case. Introducing such depth-dependent a into a solar structure model employing a non-local 
mixing length formalism results in overshooting which is approximately proportional to a at 
the base of the CZ. Although overshooting is reduced due to the reduced a, a discrepancy with 
helioseismology remains due to the steep transition to the radiative temperature gradient. 

In comparison to the mixing length models the transition at the base of the CZ is much 
gentler in the 3D models. It was suggested recently (Rempel 2004 ) that this discrepancy is due 
to the significantly larger (up to seven orders of magnitude) input energy flux in the 3D models 
in comparison to the Sun and solar models, and that the 3D calculations should be able to 
approach the mixing length regime if the input energy flux is decreased by a moderate amount. 
We present results from local convection calculations which support this conjecture. 
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1. Introduction 

Convection poses a difficult problem for stellar structure modelling. One-dimensional 
stellar structure models require a parameterization of convection in order to be able to 
yield the thermal stratification within the CZ. Since convection in the stellar envelopes 
is in general highly efficient, the stratification is close to adiabatic in much of the CZ 
and thus a detailed description of convection is not needed there. The most often used 
way to parameterize convection is to use the mixing length concept in a form introduced 
by Vitense (1953) see also Bohm-Vitense (1958). which considers convective elements 



to lose their identity after rising (or descending) the so-called mixing length which is 
proportional to the local pressure scale height, i.e. 

I = aH p . (1.1) 

Using this basic assumption, it is possible to derive equations that relate the velocity and 
temperature fluctuations to the mean stratification, and thus, to compute the convective 
energy flux (see, e.g. Chapter 6 of Stix[5002). 

Among other conceptual problems, the mixing length formalism neglects the effects of 
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Figure 1. Coriolis number in the solar convection zone according to 11.211 . 

rotation although convection can be significantly influenced by it in regions of the CZ 
where the turnover time is longer than the rotation period. In a recent study, Kapyla et 
al. I|200 5) estimated the Coriolis number, which is the inverse of the Rossby number, in 
the solar CZ from 

Co = RcT 1 = 2Q q t = 2Q Q aH p /u , (1.2) 

where a = 1.66, FIq = 2.6 ■ 10 -6 s _1 , and u the convective velocity obtained using a local 
mixing length model. This relation gives values of the order of 10 -3 near the solar surface 
and of the order of ten or larger near the base of the CZ (see figure QJ. We compute the 
mixing length parameters relating the velocity and temperature fluctuations to the mean 
stratification from local 3D convection calculations in the range Co = ... 10, which 
coincides with the range expected in the Sun. Thus it is possible to probe the influence 
of rotation on the mixing length relations and take the effect implicitly into account in a 
solar model employing a non-local formulation of the mixing length concept in order to 
study overshooting below the CZ. 

There is a striking difference between the almost adiabatic overshooting with a very 
sharp transition to the radiative gradient seen in non-local mixing length models in com- 
parison to the much more subadiabatic and smoother overshooting seen in 3D convection 
models. In a recent paper Rempel l|2004|) suggested that this discrepancy arises due to 
the fact that the two models are simply working in different parameter regimes in the 
sense that the input energy flux in the 3D models is usually up to 10 7 times larger in 
comparison to the non-local mixing length models and the Sun. We have performed 3D 
numerical calculations in which we decrease the input flux by two orders of magnitude 
in order to study this effect in more detail. 

The remainder of the paper is organised as follows: in §[21 a brief description of the 
model is given and in S I3.1I and § 13.21 the results concerning the effects of rotation and 
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input energy flux on convective overshooting are presented. Finally, the main results of 
the study and remaining problems are summarised in §0] 

2. Numerical model 

We use the same model as that described in Kapyla et al. ( 20051 l2T)()fi[l . The computa- 
tional domain is a rectangular box situated at a latitude O, in which case the rotation 
vector is represented by fl = f2o(cosO,0, — sin0). In the present study the calculations 
with rotation are performed at the south pole, i.e. = —90°. The fluid obeys the ideal 
gas law and radiation is taken into account only via the diffusion approximation. 

In contrast to many earlier studies (e.g. Brummell et al. I2002|l we do not use a piece- 
wise polytropic stratification which implies that the thermal conductivity behaves like 
a step function. Instead, we use a smoothly varying stratification where the logarithmic 
temperature gradient is computed from 

V = V 3 + i{tanh[4(z m - z)] + 1}AV , (2.1) 

where V3 = 0.15 is the gradient at the bottom, AV = Vcz — V3 the difference between 
the gradient in the unstable layer and the applied gradient, and z m the inflection point 
of the tanh- function, calculated so that V = V a d in the initial state at the base of the 
convectively unstable region at z/d = 1 . 

In order to regulate the input energy flux we split the heat conduction term, d t e — 
. . . + r conc i, in the internal energy equation into two parts 

Tcond = V • [/stV(e - e) + K h Ve] , (2.2) 

where the first term acts only on the fluctuations and the latter only on the mean, 
i.e. horizontally averaged stratification, and where e = cyT. Thus K t and Kh can be 
considered as the turbulent and radiative conductivities, which satisfy K t ^> Kh in real 
stars. We define the conductivies as 

«t = IPXo , (2.3) 
(7 - 1) Fb 

K h = = , (2.4) 

where xo is the reference value of the thermal diffusivity, computed from Pr = v/xo, 
where Pr = 0.4 is the Prandtl number and v the kinematic viscosity. Ft, is the input energy 
flux, g the constant gravitational acceleration, and V the mean logarithmic temperature 
gradient given by (|2.1|l (for more details, see Kapyla et al. 2006). 



3. Results 

3.1. Effects of rotation on mixing length relations 

Using the basic assumption of the mixing length concept, Ql.lfl . it is possible to derive 
equations that relate the velocity and temperature fluctuations to the mean thermal 
stratification 

^^f^(V-V ad ), (3.1) 

r 7I = ^(v-v ad )r, (3.2) 

where the bars denote horizontal averaging, primes the fluctuation, and V a d = (7— l)/7 = 
0.4, where 7 = cp/cy = 5/3. Furthermore, adiabatic variation within the convective 
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elements is assumed. Figure shows the results for calculations made at the south pole 
with approximate Coriolis numbers (no rotation), 1, 4, and 10. It is clear that increasing 
rotation reduces the convective efficiency and thus the mixing length parameters. If the 
rotation dependence is intepreted as a depth dependence in the solar CZ, the reduction 
of the mixing length a can be taken into account in solar models as an implicit way of 
incorporating the effects of rotation on convection. The main effect of reduced a near the 
base of the CZ is the reduction of the overshooting when a non-local version of the mixing 
length model is used (Kapyla et al. 20Q2J- The overshooting depth is approximately 
proportional to the mixing length at the base of the CZ. 

3.2. Effects of input energy flux 

The non-local mixing length models tend to produce a quasi-adiabatic overshoot region 
with a very sharp transition to the radiative gradient below the CZ (see e.g. figure 9 of 
Kapyla et al. 2005). Although rotational effects may be able to alleviate the situation, the 
sharp transition should still show up in the helioseismic inversions of the solar internal 
structure. On the other hand, numerical 3D calculations always tend to produce over- 
shooting with a much gentler transition to the radiative gradient (e.g. Brummell et al. 
2002). The main difference between these two approaches is that whereas in the mixing 
lenght models the input energy flux is the solar flux, i.e. / = Fq/ pc^ ~ 10 -11 in the deep 
layers of the CZ, the 3D models need a much higher flux (up to 10 7 times) in order to 
bring the thermal relaxation time closer to the dynamical time scale. 

Recently, Rempel f2004| suggested that if the input energy flux in 3D calculations 
was reduced by a moderate amount, the mixing length regime could be approached. Our 
results (see figure^ support this conjecture. When the input energy flux is reduced by a 
factor of 10 2 it is seen that the overshooting depth decreases as the average velocities are 
reduced, and that the transition to the radiative gradient becomes significantly steeper. 
If these results are taken at face value it would seem difficult to avoid the quasi-adiabatic 
overshoot region with steep transition at the base of the solar CZ if extrapolated to the 
solar regime which is still five orders of magnitude away in terms of the input energy flux. 
One must, however, bear in mind that in the present models the convectively unstable 
region spans only little over two pressure scale heights so the effects of compressibility 
are likely to be weak in comparison to the Sun, affecting the filling factor of downflows. 
The filling plays a crucial role in the overshoot model of Rempel (2004), with low values 
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Figure 3. Superadiabatic temperature gradient S = V — V a( j as a function of the input energy 
flux Fb. The input flux is given in units of po(gd) 3 ^ 2 , see Kapyla et al. 120U6H for the details. 
In the present models pc 3 w po(gd) 3 ^ 2 at the base of the convectively unstable region leading 
to dimensionless flux of / = F^/pc 3 « 0{F\,). The thin dotted curve shows the temperature 
gradient if the total flux would be transported by radiative diffusion. 

(w 10~ 5 ) being able to produce overshooting with smooth transition also for the solar 
energy flux. 



4. Conclusions 

Three dimensional local convection calculations were used to probe the effects of rota- 
tion on mixing length coefficients relating the temperature and velocity fluctuations to 
the mean thermal stratification. It was found that when the rotational influence on the 
flow is comparable to that expected in the deep layers of the solar CZ, the mixing length 
parameters are reduced by a factor of three to five. If a depth-dependent mixing length 
a is introduced into a solar model, the overshooting at the base of the CZ is reduced 
approximately in proportion to the reduction of a. Although the depth of the solar CZ 
can be correctly reproduced in this way, the steep transition to the radiative gradient 
should still be visible in the helioseismic inversions. 

The overshooting in 3D convection calculations is much more subadiabatic with a 
smooth transition, which is due to the much higher input energy flux used to meet the 
time step constraints. In the present study we show that decreasing the input flux in 3D 
calculations leads to more adiabatic overshooting and sharper transition at the base of 
the CZ. Although this result seems to suggest that the 3D calculations will approach the 
mixing length regime when the flux is reduced enough, one must bear in mind that in the 
present models the stratification is rather weak in comparison to the Sun. Thus the effects 
of compressibility are likely to be underestimated, and lead to too large filling factor for 
the downflows. Furthermore, the spatial size of the downflow plumes is restricted by the 
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grid size, and it is probable that the filling factor of downflows further decreases when 
the resolution is increased. 

To summarize, we stress the point that rotation should be taken into account in models 
of convective overshooting since it can exert considerable influence on convection already 
in slowly rotating stars such as the Sun. Furthermore, high resolution numerical studies 
of deep convection are needed in order to study whether the convective overshooting is 
due to very few strong downflows, producing nearly adiabatic overshoot region with a 
steep transition, or whether downflows of different strengths penetrate into the stable 
region in a larger area producing smooth overshooting required by helioseismology. 
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Discussion 

J. Christensen-Dalsgaard: Comment: Helioseismology shows that the sound speed 
gradient is likely smoother in the Sun than in models even without penetration, possibly 
requiring a subadiabatic gradient in the lower parts of the convection zone. 

R.F. Stein: Comment: when you model the base of the convection zone, it is very 
slightly subadiabatic inside the convection zone. 

P.J. Kapyla: This is indeed the case also in our 3D models. 



